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We present a tomographic approach to the study of quantum nonlocality in multipartite systems. 
Bell inequalities for tomograms belonging to a generic tomographic scheme are derived by exploiting 
tools from convex geometry. Then, possible violations of these inequalities are discussed in specific 
tomographic realizations providing some explicit examples. 


I. INTRODUCTION 

The Bell inequalities |Tj demonstrate paradigmatic difference of quantum and classical worlds. They 
were originally written for dichotomic (spin—variables O. Spin—| operators realize the Lie algebra of 
the SU(2) group. For several spin particles their spin operators form Lie algebra of the tensor product of 
the Lie algebras. Due to algebraic equivalence of the operators satisfying commutation relations of the Lie 
algebra constructed from particle spin operators and constructed from creation and annihilation operators of 
a field, one can obtain Bell inequalities also for the case of continuous variables besides discrete ones f3]. 
Beyond the specific operators involved in fhe Bell inequalifies, fheir possible violafions obviously depend 
on fhe slate under considerafion. 

For a (mulfiparlile) classical syslem wifh fluclualions, fhe syslem slafe is described by means of a joinl 
probabilily disfribufion function of random variables corresponding to fhe subsyslems. In conlrasl, for a 
(mulfiparlile) quanlum syslem Ihe slate is described by Ihe density malrix. In view of Ibis difference Ihe 
calculations of Ihe system’s slalislical properties (including correlations) are accomplished differenlly in 
classical and quanlum domains. 

Recenlly, a probability represenlalion of quanlum mechanics has been suggested |4|. This represenla- 
lion, equivalenl to all olher well known formulations of quanlum mechanics (see, e.g. |j5i|), goes back to 
quantum tomography, a technique used for quanlum slate reconslruclion ||6l. The approach makes use of a 
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set of fair probabilities, tomograms, to “replace” the notion of quantum state. It has also been understood 
fT] that for classical statistical mechanics the states with fluctuations can be described as well by tomo¬ 
grams related to standard probability distributions in classical phase-space. A comparison of classical and 
quantum tomograms can be found in Ref. 171 HI. 

Thus, in the probability representation, tomograms turned out to be a unique tool to describe both classi¬ 
cal and quantum states. As a consequence they represent a natural bed where to place inequalities marking 
the boarder-line between quantum and classical worlds. Tomograms can be either continuous or discrete 
variable functions depending on the tomographic scheme (realization). In both cases they might be directly 
used to test nonlocality. This possibility was described for symplectic tomography Q in bipartite system 
ITOl and spin tomography ifTTl still in bipartite system ifTTl . 

Here we shall derive Bell inequalities for multipartite systems in terms of tomograms belonging to a 
generic tomographic scheme. Then, we shall discuss the possibility to violate such inequalities depending 
on the tomographic realization. 

The layout of the paper is the following. In Section]^ we formalize quantum tomography in a multipar¬ 
tite setting. Then, in Section III we derive the Bell inequalities in terms of tomograms. In Section IV we 


provide some evidences of violations of such inequalities for spin—| systems as well as for field modes and 
finally draw the conclusions in Section [V| 


II. QUANTUM TOMOGRAPHY 

Here we briefly review the general quantum tomography approach for a single system, by detailing three 
relevant cases (optical llT3l . spin lITTIl and photon-number tomography HU) and then extend the formalism 
to multipartite systems. 

The basic ingredients of any tomographic scheme are a Hilbert space V. associated with space of the sys¬ 
tem under consideration and a pair of measurable sets {X, A) with measures p{x) and correspondingly. 
More precisely, the set of system states is the set S{'H) of Hermitian non-negative trace-class operators on 
H with trace 1. Usually the set X is the spectrum of an observable of the system and the set A plays the 
role of transformations. 

We use the notation V{X) for the set of probability distributions on X, i.e. the set of nonnegative 
measurable functions p : V —)• M normalized to one in the following sense f p(x) dp{x) = 1. 

Both sets S{'H) and V{X) are closed with respect to the convex combinations: if p, d G S{T-L) (resp. 
p(x), q{x) G V{X)) and a G [0,1] then 


aQ+ {1 — a)(T G S{T-L) (resp. ap{x) -h (1 — a)q{x) G V{X)). 
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Definition 1 A map T : S{'H) ^ is called tomographic map if the following three conditions are 


satisfied: 

1. for any g G S(T-L) the image T{g) ■ X x A ^ M. restricted on the set X x {A} w a probability 
density on X 


Tx{g)£V{X) VAgA, where rA(^) = r(^)Ux{A} : ^ ^ K. 


2. the map f preserves convex combinations 


T{ag + (1 — a)d) = aT{g) + (1 — a)T(cT), V^, d G S{'H), a G [0,1]. 


3. the map T is one-to-one 


T{g) = T{d) 4^ g = a. 


These conditions have simple meaning: (i) means that the tomogram T{g) of any state p is a probability 
distribution on X parameterized by the points of A, (ii) is the linearity condition, and (iii) requires that the 
tomogram of each state be unique, or, in other words, that any state can be unambiguous reconstructed from 
its tomogram. 

In the present work we deal with tomographic maps of the following form 



( 1 ) 


where U{x, A) is a family of operators on H parameterized by points (x, A) of the set X x A. In the 
examples considered below the state g can be reconstructed from its tomogram p§{x, A) according to the 
formula 


Q = 



XxA 


p{x, X)'D{x, A) dp,{x) dv{\) 


( 2 ) 


for the appropriate (x, A)-parameterized family of operators ^(x, A) on %. 


The set X is the spectrum of an observable O and the set A is a group equipped with a representation 
(in general projective) tt ; A —)• in . The operators U (x, A) have the following form 


U{x, A) = 7r(A)|x)(x|7r'^(A), 


(3) 


where |x) is an eigenstate of the observable O. For a group theoretical approach to quantum tomography 
see llTSl . See also |[T6l for a relation to grupoids. 
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A. Spin tomography 

Let us consider a system with spin j. In this case we have: T-L = X = {—j, — j + 1, ..., y — 1, j} 

and A = SO(3,M). We denote the elements of the sets X and A as s and respectively. The measure 
on X is equal to one on each element, so the corresponding integral is simply the finite sum over 2j + 1 
terms. The measure on SO(3, M) is Haar’s one. For the group SO(3, M), parameterized with Euler angles 
= {ip, Tp, 6) the measure reads 1/(11) = v{ip, ip, 9) = simp dp dp) d6 and the operator 17 of ([^ takes 
the form 


U{s,Q) = k{n)\j,s){j,s\KH^)- 


(4) 


Here the vectors \ j, s), s = —j, — j + 1,..., j — 1, j are the basis of the space (eigenvectors of the 

spin projection Sz) and the operators K{Q) are the operators of the irreducible representation of SO(3, M) 
in Their matrix elements are given by 


X cos^~^^'{p)/2) sin^'{pj/2)P^p~k^ '^^\cosp)), 


(5) 


with the Jacobi polynomials. 

Then the tomogram p(s, H) = p{s, p, pj, 9) of 0 is 

p{s,^) = {j,s\k{n)Qk^n)\j,s) 


( 6 ) 


Due to the property {j, s|iT(ll)|j, s') = (—1)'^ —s|.^(ll)|j, —s'), the tomogram does not depend on 

the angle 9, i.e. p{s, p, pj, 9) = p(s, p, p)). 

Finally, the operator D of (|^ results 

3 

V{s,Q.)= ^ (y,n|D(s,ll)|j,m)|j,n)(j,m|, 

n,m=—j 

where the matrix elements {j, n\'D{s, ll)|j, m) are given by the following expression 


(y,n|D(s,ll)|j,m) 


Svr^ 




j3=0 


h 

X ^ (j,A:|iT(ll)|j,0) 

k=-j3 


j j h 

n —m k 


I j j J3 
\ s —s k 


(V) 


in terms of Wigner 3j-symbols. 
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B. Optical tomography 


( 8 ) 


Here we have : % = L 2 (M), X = M and A = {e*®|0 G [0,27r]}. The measures on X and A are 
Lebegue’s ones. The operator corresponding to Eq.Q reads 

ij{x,e) = R{e)\x){x\R\e), 

where R{6) is the rotation operator 

R{e) = exp , 

acting on and the canonical position x and momentum p operators as 


. X \ I cos 9 — sin 9 

R{9) I ^ R\9) = 

pi \ sin0 COS0 


X 

p 


In other words, U{X,9) of ([^ is the projector on the rotated eigenvector \X) of the position operator x. 
The tomogram p{X, 9) of 0 is the diagonal matrix element 

p{X,9) = {X\R{9)eR\9)\X). (9) 


dr. 


Furthermore, the operator 2? of (M results 


1 


/' 


T){X,9) = — I |r| exp( —ir(X — cos— sin 


C. Photon-Number tomography 


Here we have: R = L 2 (M), X = Z+ = {0,1,...} and A = C. We denote the elements of the sets X 
and C as n and a respectively. The measure on X is equal to one on each element and the measure on C is 
(l/Trjd^a, where d?a = dReadIma is the Lebegue’s measure on the real plane. Here, the operator IJ is 
the projector onto the displaced Fock state 

U{n,a) = D{a)\n){n\D^a), (10) 

with 


D{a) = exp 

From ([T]l the tomogram p(n, a) reads 


a — a ^ a + a X. 

X — ^^—p 


v/2 


v/2 


p{n,a) = {n\D{a) gD\a)\n). 

Furthermore, the operator P of 0 becomes in this case 

+ 0O 

V{n, a) = 4(-l)" ^ {-irD{a)\m){m\D^{a). 

m=0 


( 11 ) 
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D. Tomography for multi-partite systems 

The generalization for multi-partite systems is straightforward. 

Definition 2 Consider a n-partite system with the state space and n tomographic schemes, one for 
each part with sets {X^, A^) and operators tJk{xk, Xk) and Vkixk, Xk), k = 1,..., n. The tomographic 
scheme for the whole system is then constructed as the direct product of these schemes, by using 

n n 

X=llXk, 

k=l k=l 

n n 

U{x,X) = (^Ukixk,Xk), f){x,X) = (^Dkixk,Xk), ( 12 ) 

k=l k=l 

where x = (xi,..., Xn), A = (Ai,..., Xn) and the measures p{x), u{\) on X, A are direct products of 
pi{xi ),..., Pn{xn) and z/i(Ai),..., r'n(An) respectively. The tomogram p{x, A) of a state g (generalizing 
Q)is 

p{x,\) = TX:(^QU{x,\)y (13) 

For any X € Ait is a probability distribution on X, thus p(x, A) dp{x) = 1. 

Remark. From the definition ( [l^ of the operator U{x, A) it immediately follows that the tomogram p(a;, A) 
of a factorized state 


g = §10 ...iSi Qn 


(14) 


is also factorized, i.e. 


p{x,X) =pi{xi,Xi) . . .pn{Xn,Xn), 


(15) 


where pk{xk, Xk) is the tomogram of the state gk- More generally, the tomogram of a separable state 


+ 00 

M) 

g = > aig\' 


X 

i=0 


fi) 


-1-00 


ai ^ 0, ^ ^ Oi — 1 


(16) 


j=0 


is also separable in the following sense 


+ 0O 


p{x,X) = '^aiP^^\xi,Xi) . . .pW(Xn, An), 


(17) 


i=0 


where (xk, Xk) is the tomogram of the state 
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III. BELL INEQUALITIES FOR TOMOGRAMS 


Let us consider a n-partite system in tomographic representation, whith each subsystem supplied by a 
tomographic map 7fc, k = 1,... ,n. The tomogram p{x, A) of a state p is a function of 2n arguments and 
with respect to one half of them it is a probability distribution. We will show that in general it cannot be 
considered as a classical joint probability. 


Definition 3 For any k = 1,... ,n let and Zk be two measurable sets such that 


Xk = Yk\JZk, = 

and for any Xk 6 let Ak{Xk) be a dichotomic random variable on X = such that 

P(^fc(Afc) =-1) = fz^TT(^gUk{xk,xf dpkixk)- (18) 


Symbolically the variables Ak{\k) can be written as 


Ak{Xk) 


1 A Xk e Yk, 
-1 if Xk G Zk 


(19) 


in the coordinate system deformed by the operator Uk{xk, Xk)- The joint probability distribution of the 
random variables yli(Ai),..., An{Xn), namely 


Pei,...,en (^1) • • • ) Xn) — P(^l (Al) — Si,..., A^ii^Xn) — £^n) > 


where Sk = ±1, is given by 


with 



p{x, A) dx, 


Wk 


Yk if Ek = 1, 
Zk if Ek 1* 


The correlation function of Ai(Ai),..., An{Xn) results 


^(Ai,...,A„) = (yli(Ai)...A„(A„)) 

~ ^ ^ P£l,...,£„ (^1) • • • )-^n)^! • • • Sn- 

ei,...,£„=±l 


( 20 ) 


( 21 ) 



Definition 4 Let us fix two parameters and for k = 1,..., n and denote 


E{ju ..., jn) = E{X^i^\ ..., jk = 1,2. 


( 22 ) 


Since each index jk can take 2 values independently on all the other indices, there are 2^ correlation 


functions (221. Then we define by 

the vector of these correlation functions with some order of multi-indices {ji,..., jn)- 


(23) 


It is convenient to enumerate the functions E{ji,..., jn)- For this purpose we use the binary base with 
“digits” 1 and 2 instead of 0 and 1. This means that we use the following one-to-one correspondence 

{1,..., 2 } 9 y (j'l, ... ,jn), jk = 1) 2, 


where j and (y'l,..., jn) are related to each other according to 

j = {jl ~ 1)2"' ^ {jn — 1)2 -|- jn- 

By virtue of such an ordering, the vector e (231 can be written as 

e= (i^(l),...,i?(2")) = (i^(l,...,l),...,S(2,...,2)) 


(24) 


(25) 


What region Qn C fills the vector e of (251? Due to the fact that each observable has only two 


outcomes ±1 it follows that each correlation function (221 is bounded by one by absolute value and, so, the 
set Lin is a subset of 2”^-dimensional cube [—1,1] . 

Suppose that it is possible to model the result of the measurement by a random variable, Ak{jk), which 
can take two values ±1. We assume that these random variables can be arbitrary correlated. 

Definition 5 Let us define by 

p(zi(l), ...,Zn(2)) = p(Ali(l) = zi(l),...,A„(2) = z„(2)), (26) 

the joint probability distribution for random variables Ak{jk), with ik{jk) = 3=1. Since each index ikijk) 


can take independently 2 values we have 2^^ numbers (261 which completely describe statistical character¬ 
istics of the random variables under consideration. WF enumerate them with a single number i = 1,., 2^” 
using the same rule as for the correlation functions E{j), namely 


{l,...,2^^}3i^{n{l),...,in{2)), 
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where i and (ii (1), • ■ •, *n(2)) are related to each other according to 

i = (zi(l) - 1)22-1 ^ + i^(2). 

Enumerated in such a way the probabilities \26\form a 2^^-dimensional vector 

/ \ ,,. r)2n 

P= {pi,...,P22n) GW . 


(27) 


(28) 


The point (281 lies in the standard simplex 


522n_l — < (pi, . . . ,P22n) 


'^Pi = I,Pi ^ 0 I C 


(29) 


i=l 


What region C flPs the vector e (25) when the point p (28) runs over the simplex 522n_ 


I? To answer this question we are going to explicitly relate e and p assuming the former expressed like 
classical joint probabilities. Then, the correlation function E{j) is intended as a simple linear combination 
of Pi with proper coefficients. Looking at ([IT]) we consider such coefficients, given by the product 


= hijl) ■ ■■inUn), 


(30) 


where jj- and ikijk) are “digits” of the numbers j and i in the binary representations (24) and (27). The 
numbers S{j, i), j = 1,..., 2”, z = 1,..., 2^” form a 2"^ x 2^”^ matrix Sn and the relation between e and p 
can then be written as 


e = EnP- (31) 

We see that the region is the image of the standard simplex S 22 n_i 

= Sn{S22n_i), (32) 

where we do not distinguish the linear map Sn ■ —)• and its matrix Sn in the standard bases of 

M2^"andM2". 

Thus, we have reduced the problem of finding Bell inequalifies fo find fhe sef Qn- It means fhaf fhe 
problem of finding Bell inequalifies boils down fo a sfandard problem of convex geomefry, referred fo 
as convex hull problem: given poinfs Cj find fheir convex hull, or facefs of maximal dimension of fhe 
corresponding polyfope (for nofions of convex geomefry see, e.g. ifTTll i. 

Now we will gef fhe Bell inequalifies explicifly. Nofe thaf permufafions of fhe columns of the matrix Sn 
do not change their convex hull and that they correspond to permutations of the components of the vector 
p or different orderings of the probabilities ( [28] ), so one can safely permute columns of Sn without altering 
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Theorem 1 The set Tin specified by the (Bell) inequalities for the vector of the correlation functions 


(e, H 2 nc) ^ 2"', Vc = (±1,..., ±1). 


The matrix i/ 2 ” i^ the Hadamard matrix recurrently defined as 


(33) 


i/2” — ^2 


H2, H2= 


1 1 


1 -1 


Proof. The key fact in deriving the Bell inequality (33 1 is that the matrix Sn can be written in the following 
block form 


(34) 


Tn — ( i/2^ —i/2” • • • i/2” —i/2” 


after appropriate arrangement of its columns. One can rewrite the r.h.s. of (34i as the product of two 
matrices 


S-n — i/2” ^ ii/2” —£'2” • • • i?2” —£'2” J — .^2”-4ri; 

which means that the linear map £n ■ —)■ can be decomposed into two maps 


Sn — i/ 2 ” C 




■n on 


According to this decomposition ((3T]) reads 


e = // 2 ”q, 


where the vector q = AnP G is explicitly given by the following expression 


/ Pi — P2”+l + ■ ■ ■ — P(2”-l)2”+l ^ 


q = 


y P2” — P2-2” + . . . — P22” / 

Define the following convex polytope Oat C 


Om = {x ^ M^K®, c) ^ 1, Vc = (±1,..., ±1)}. 


(35) 


(36) 


(37) 


As one can easily see the image of the standard simplex 522n_i is exactly the polytope 02”, that is 
An{S 22 n_i) = C> 2 ”- From this fact we have 


Tin — i/2”(672”)- 


(38) 
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Now the Bell inequalities can be straightforwardly obtained from this relation. Just notice that a non¬ 
degenerate linear map / : —>• with the matrix F maps a half-space f) = {s G a) ^ b} to 

the half-space /(ft) = {y G a) ^6}. Taking into account the following representation of 

the polytope ( [3^ 

^2" = n {9l(9,c)^l}, (39) 

c=(±l,...,±l) 

the symmetry of the Hadamard matrix H 2 " and the formula for its inverse H^n = we get the 

explicit form of the set i.e. 


Pi {e|(e,f?2-c) ^2’^}. 

c=(±l,...,±l) 


Hence, the Bell inequalities (331. 


(40) 


Remark. Explicitly ( [3^ can be written as 

2 


■ ■ ■ , jn) 

ilv,in = l 




(41) 


where the coefficients are connected with the vector c by the following relation 


Jn = ^ c(ei,..., £n)e’i ^ (42) 

The number c(ei,..., £n) here is the z-th component of the vector c, where the binary representation of i is 
z = (ei... £ 11)2 with digits 4-1 and —1 instead of 0 and 1. 

One can easily see that there are 2”+^ inequalities of the form ±£^(ji,...,/„) ^ 1. They correspond 
to the functions c(ei,..., e„) that are columns of either or —Tf 2 " and they are referred to as trivial 
inequalities. 

Finally, notice that the well known CHSH inequality ifTFI is a particular instance of ( |4T| ) corresponding 
to n = 2 and 


c(—1, —1) = — 1 , c(— 1 , -fl) = c(+l, —1) = c(+l, -fl) = -|-1- 

Theorem 2 Any separable state satisfies ( |33p with correlation functions ( |22p . 

Proof. Let us start with a factorized state ( [T4| whose tomogram ( [T5] ) is also factorized. Due to this the 
random variables ^i(Ai),..., An{Xn) are independent and the correlation function E{Xi ,..., An) reads 


E(^Xi, ■ ■ ■ , An) — 9l('^l) ■ ■ • QniXn) 


(43) 






where 


pi’lK^k) = PiAki^k) = £k), = ± 1 - 


Due to the fact that 


Pi^\^k)+P-l{>>k) = yk = l,...,n VAfc G Afc, 


it is clear that —1 ^ qk{^k) ^ 1- The left hand side of the inequality (331 is a linear function of any 
Qki^k'^^) where all the , qn(An"^), jk = 1, 2 are considered as independent variables. A linear 

function defined on the convex set [—1,1] takes its maximum on a boundary point, ±1 in this case, and so. 


the left hand side of (331 is maximal if ~ Jfc = 1) 2, fc = 1,..., re. In such a case the vector e 


of correlation functions is a column of either 1^2" or —Just note that due to (431 the vector e reads 


e = 




where qk{jk) = qk{^k'‘^)- ^ = Q is the i-th column of H 2 ^ or —Tf 2 ", then 


(e,if2"c) = ±(cj,Tf2"c) = ±(i/2"Cj,c) = ±(2"'ei,c) = ± 2 ^ ^ 2 ^ 


(44) 


Here we used the orthogonality of the columns of H 2 r>.: H 2 «Ci = 2”ej, where all the coordinates of e* are 
zero except the z-th which is one. Hence, we have proved that all factorized states satisfy ( |^ . 

Let us now consider a general separable state ( [T^ . Since the correlation function E{\i ,..., A^) is a 
linear function of the state, the vector e is a linear combination of the vectors corresponding to the 
states (g)... (g) Qn \ i.e. 


+OD 

e = Oje 

n=0 


(i) 


As we have already shown each vector satisfies all the Bell inequalities (33 1 or lies in the convex set 
Once all the vectors are in so is their convex combination e. This means that any separable 


state satisfies all fhe inequalities (331. ■ 
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IV. QUANTUM VIOLATIONS 


The Bell inequalities are of interest not because they are always valid but because they can be violated. 
One can ask if there was a mistake in the proof of theorem [T] The problem relies on the underlying hypoth¬ 
esis of locality when relating e with p in pT] ). In doing so we have implicitly assumed ( [2T] ) as a classical 
joint probability, which is not generally true at quantum level. 

We follow Mermin lfT9l to derive the only Bell inequality whose maximal quantum violation is the 
largest among all others. 

For an odd number n of systems let us consider the following random variable 


Mn = Im 


.k=l 


(45) 


Since each can take only values ±1, each term in this product is equal to y/2 by absolute value. Further¬ 
more, since n is odd the whole product has the phase that is an integer multiplier of vr/d. Asa consequence 
we have 




(46) 


Explicitly this inequality reads 




< 2 (’"- i )/2 


(47) 


where the sum here runs over the set of multi-indices (ji,..., jn) which contain an odd number of 2 

J = {(ji, • • •,in) = 2}| = 2/ -h l| 

and 


• ■ ,jn) = I, \{k\jk = 2}| =21 + 1. 


Multiplied by the inequality (471 takes the form (41 1 and it is easy to show that it is a Bell inequality. 


i.e. there is a vector c that gives the coefficients of (47 1 (multiplied by 2("'+^(/^) according to (42 1 . 


We now consider an even number n. Let us denote the expression (451 as M„(l, 2) and the similar 


expression with the variables Afc(l) and Ak{2) swapped as M„(2,1). Consider the following combination 


Mn = 2)(A„(1) + An{2)) + M„_i(2, 1)(A„(1) - An{2)). 

Since 2) is equal to ±2"'/^“^ and An{j) = ±1, we will have 


(48) 


(49) 
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Using the explicit form (|47|) for the odd number n — 1 one can write (49 1 as 






where 


Kh, ■ • • , jn) 


1 if jn = 2 and |{A;|jfc = 2}| is nonzero and even 
0 otherwise 


(50) 


One can see that it is a Bell inequality and multiplied by 2*^/^ it takes the form (411. Furthermore, forn = 2 
Eq.(50l exactly reduces to the CHSH inequality |[T8]| . 

Let us now see how the inequalities ( [47] ) or (50) can be violated in different tomographic realizations 
starting from the following entangled state 


= 72(10)+ (51) 

where 0 = (0,..., 0) and 1 = (1,..., 1). Below, for the sake of simplicity, the focus will mainly be to 
n = 2,3. 


A. Spin tomography 


Using the notation |0) = [ — 2 )^ |1) = I + 2 ) for tho spin projection along z, the state (511 becomes 


1 


1 '^') — • • • > “ 5 ) + l + 2 > • • • ’ + 2 )) 


whose tomogram, referring to Eq.(0, reads as 

1 


p(si, . . . , Sri) • • • ) ffn) — 2 


j=i i=i 


For n = 2 we immediately get 


F (+2 > +2 ’ ^ 1 ’ ^ 2 ) 


p( + |, -|,Ui,U2) 


= ^(1 + cos7i 008 72 + sin7i sin 72 cos((^i + (^ 2 )), 
= “ cos7i cos 72 — sin7i sin 72 cos((^i + (^ 2 )), 


(52) 
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and the correlation function (211 becomes 


£'(^1,122) = cos V'l cos 1^2 + sin'll sin■02 cos((/?i + ^ 2 )■ 


The Bell inequality (501 reads in this case 


^ 2 , 


for all 12^'^ = j,k = 1,2. The maximum of the l.h.s. of (541 with (53) 

attained by taking e.g. (the angles 6 do not matter here) 

= (<^i,-7r/8,0), 12^^^ = (-v?i,7r/8,0), 

f2f^ = ((^i,37r/8,0), 12^^^ = (-(^i,-37r/8,0). 


In the case of n = 3, from (521, we have (not to overload the notation we omit the f2’s) 


p( + 2 , + 5 , + 5 ) = |[1 + cos V'l cos-02 + cos'll cost's + cos-02 COS'03 

- sin'01 sin 0)2 sin03 cos((^i + (p2 + y^a)], 

p{ + \, +2’ “ 5 ) ~ + cos 01 cos 02 — cos 01 COS 03 — COS 02 COS 03 

+ sin01 sin02 sin03 cos(y)i +ip2 + p;i)\, 
p{ + \, “5) +5) = |[1 “ cos 01 cos 02 + cos 01 COS 03 — COS 02 COS 03 
+ sin 01 sin 02 sin 03 cos((^i + ip2 + ¥^3)], 
p{ + \, —5, —5) = |[1 — cos 01 cos 02 — cos 01 COS 03 + COS 02 COS 03 

- sin01 sin 02 sin 03 cos(y)i +ip2 + y^s)], 

p(— 2, +2’ “l"2) ~ “ cos 01 cos 02 — cos 01 cos 03 + cos 02 cos 03 

+ sin 01 sin 02 sin 03 cos(yPi + ip2 + ¥^3)], 
p(—2, +2’ “s) ~ “ cos 01 cos 02 + cos 01 cos 03 — cos 02 cos 03 

- sin 01 sin 02 sin 03 cos(vpi + ¥)2 + ¥^3)], 

p(—2, “5) +^) = |[1 + cos 01 cos 02 — cos 01 cos 03 — cos 02 COS 03 

- sin 01 sin 02 sin 03 cos(vpi + (/92 + ¥^3)], 

p(—2, “2’“5) + cos 01 cos 02 + cos 01 cos 03 + cos 02 cos 03 

+ sin 01 sin 02 sin 03 cos((/?i + ¥^2 + V^a)] ■ 

Thanks to these tomograms the correlation function (|2T]) results 


(53) 

(54) 

is 2\/2 and is 


£(ffi, 02,03) 


sin 01 sin 02 sin03cos(¥9i + p 2 + P‘i)- 


(55) 
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Finally, the Bell inequality ( [47] ) in this case reads 


< 2. (56) 


Using (55 I the maximum violation occurs when the l.h.s equals 4. This value can be attained by taking e.g. 
(again the angles 9 do not matter here) 

= 7r/2, = 57r/6, 

= 7r/2, = iff'' = = 7r/3. 


B. Optical tomography 

The tomogram of the state ( [5T] ) accordingly to ([^ is given by 

p{X, 6) = 1 + 2^ [](X2) + 2 ("+ 2)/2 ]^(Xi) cos(0i + ... + 0, 


exp 


i=l i=l 

where X = {Xi ,..., X^) and 0 = (0i,..., On)- 
We take the sets and of Definition]^ to be 

Yfc = [a:, +CX)), Zj, = {-oo,x). 

For such sets and tomogram (57 1 , the correlation function ( |2T] ) results 

E{e) = 2^-1 [a^(x) + ay(x)] + 2"b^(x) cos(0i + ... + 9n), 

where 

1 1 1 _ 2 1 _ 2 
ao(x) = --erf(x), ai(x) = --erf(x) + -^xe ^ , bo(x) = -j=e ^ . 
2 2 V27r 


i=l 


(57) 


(58) 


We have now to insert ( |58] ) into ( |47] ) or (501 to get an explicit version fo the Bell inequality. In doing so 
we use a Lemma, reported in|^ showing that the maximal value of 

Yj “ii.-dn cos{e^Y + • • • + 

does not exceed and this value is attained with coefficients of or ( |50| . It then follows that 

the maximal value /n(x) of the l.h.s. of ( |4^ and of (501 is 

2"|a[((x) + a"(x)| + 2’"+("+i)/2|b"(x)|, n odd 
2"'|aQ (x) + a”(x)| + 2"’+"'/^ |bQ(x)|, n even 

Figure illustrates the function /n(x) for n = 2,3. A tiny violation of Bell inequality only occurs for 
n = 3. 


fn{x) = 


(59) 
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X 


FIG. 1; Function /„ of Eq. (59 1 versus x for n = 2 (dashed line) and n = 3 (solid line). 


C. Photon-Number tomography 


Considering the state (|5T]) its number tomogram ([TT]) can be computed as 


p{mi,...mn,ai,...an) = OILi ""'LTl ^ lOILi «* + nr=i("ii “ l«*P)| • 

We further choose the sets of Definitionas Zi = ... = Zn = {0}, Yi = ... = = {1, 2,3,...}. 

The corresponding correlation function (|2T]) forn = 2 is 


(60) 


77 ( 01 , 02 ) = e [2 + 4517 ( 0102 ) + 2 |aip|o 2 '^ 

- (1 + 102^ - (1 + |«i|2) 


(61) 


Furthermore, the Bell inequality for the number tomogram with n = 2 is from (501 




^2, 


for all o^'^^ o^'^^ G C, j = 1, 2. Figure [^illustrates that this inequality can be violated using (|61|. 


Analogously, from (60 1 it follows that the correlation function (|2T]) for n = 3 


IS 


77(oi, 02 , 03 ) = e [—4 + 8517 ( 010203 ) — 4|oip|o2p|o3| 


+2(el“il^+el“2l^ + el“3l 

|2U„|2J«3|^ 


+ 2|o2|^|o3pel“^l^ 
|2UJ2 


+2|aino2f + 2|oino3fe'' 

- (1 + |oip) el"2p+l“3p - (1 + |a2p) 

_ _l_ |o3p) + el"il^+l«2p+|a3|^ 


(62) 


(63) 
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0.0 0.5 1.0 1.5 2.0 2.5 3.0 


a 


FIG. 2: The left hand side of \62\ as a function of (dashed line); the other parameters are given by = 0.165, 


= —0.165, = —0.559. The left hand side of (64 1 as a function of (dashed line); the other parameters 


are given by = 0, = 5.936, q;[^^4.767, = 4. 


This time the Bell inequality for the number tomogram reads from ^T\ 


hi{p'^ ,Oi2 , Ct3 ) + ,02 i(^2 ) ~T ,0.2 ,02 ) J^[Oi ,02 ,02 ^ 


^ 2 . 


(64) 


This inequality, by numerical checking, results never violated with (631 and an example of the behavior of 
the l.h.s. is shown in figure]^ 

By also choosing Zi = ... = = {0,..., m}, Yi = ... = Yn = {m + 1, m + 2,...}, with m > 0, 

neither (|62|) nor (641 will result (by numerical checking) ever violated by using ([^ and ([63]) respectively. 


V. CONCLUDING REMARKS 


As we have seen from the previous examples the use of finite (namely 2*^) number of tomograms within 
a tomographic realization may lead to the evidence of nonlocality. Actually, it results that finite dimensional 
systems by means of spin tomograms allow for the best evidence of nonlocality. In contrast, violations of 
Bell inequalities seem much harder to uncover in infinite dimensional systems where V. = L 2 (M). Given 
that we have considered in both cases the same (entangled) state ( [ST] ), this difference, according to Ref. EOll . 
must be ascribed to the diversity of observables employed (from which the tomograms stem). However, we 
argue that also the way the spectrum of an observable is binned could play a role. As matter of fact the 


choices made in Sections IV B and IV C for Y^ and do not exhaust all possibilities of these measurable 


sets. Unfortunately looking at Bell inequalities violations using optical tomograms (resp. photon number 
tomograms) by scanning the possible sets Y^ and Zk appears a daunting task. 
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All in all the advantage of the tomographic approach is to allow to to find the large violations of Bell 
inequalities typical of spin systems also in infinite dimensional systems. In fact, introducing in L 2 (M)®’^ 
the following local pseudo-spin operators |[^ 

H-OO 

= ^{\‘^nk){2nk + l\ + \2nk + l){2nk\^, 

nfe=0 

+ O0 

= -z ^ (^\2nk){2nk + 1| - \2nk -f l)(2nfc|^, 

nfe=0 
H-OO 

= y;(-innt)(n*|, 

nfe=0 

where \nk) are Fock states of the kth subsystem, we can derive the tomograms of the spin tomography 
realized with the above operators from those of any other tomographic scheme (see e.g. Il22]l l. The price 
one ought to pay in such a case is the completeness of the set of starting tomograms, (i.e. a number of 
tomograms much greater than 2”). 


c{k) 

Dy 

g(k) 
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Appendix A 


Lemma 1 For any coefficients aj = (j = {ji, ■ ■ ■ ,jn)) of (42) and for any angles 6^^'^ 

(k = 1,... ,n) we have 

2 


1 


Qj cos -I- ... -I- 
j=t 




The equality is attained with coefficients from (47), (50). 


Proof. To estimate the l.h.s. of ( |A1| ) note that 

2 

aj cos (^9^^^ -|- ... -h 9^"^ 
i=i 




E' 

7=1 




so we need to estimate the last sum. To this end we use ( [42] ) obtaining 
2 




7=1 


ei,...,e:n=il k=l 


ie 


( 1 ) 


ie 


( 2 ) 


(Al) 


(A2) 


(A3) 
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Next we define 


3 ( 1 ) _ ^( 2 ) 




a _ "k "k _ ^k ' 

“k — -X-, ^k — 


(A4) 


SO that the r.h.s. of (|A3|) simplifies to 


2ng*(^l+...+c,„) ^ . . . , en) n ^kiSk), 

S'!,... ,£ti— il 


k=l 


where afc(+l) = cos 9k and Ofc(—1) = i sin 9k- Taking into account that we use absolute value in (A2 1 and 
divide by 2” in (|A1|) we have to prove the following inequality 


c{ei,...,en)Y\akisk) 

ei,...,£n=±l 


k=l 


< 2("-i)/2 


(A5) 


for any ±l-valued function c(ei,..., Sn)- We employ the induction method. For n = 1 we simply have 


c(+l) cos 01 + c(—l)isin0i 


= 1 = 


Then, we can write the sum in (A51 as 


Y c(ei,... ,en) afc(efc) = A„_i cos0„ + sin0„, 

£l,...,£n=±l 


k=l 


n —1 


Y c(ei,... ,en_i,+ 1 ) afc(efc) COS0,, 

£^l5-.*,£^n— 1 —il 


k=l 


n -1 


+i Y c(ei,...,en-i,-l) afc(efc)sm0„ 

£l,...,£n_l=±l 


k=l 


where, according to the induction assumption, we have 

\A 1 1 I R 1 1 < 

Kin—1 ) \-Dn—l\ ^ ^ 


(A6) 


The sum in (A5 I can be estimated in the following way 


Y c{ei,...,en)Y\ak{ek) 

£1 ,. ,£ti — 


k=l 


= An-i cos 9n + iBn-i sin 9n 

^ V|A„_i|2 + |R„_i|2 ^ 


Now we show that with the coefficients of (47 1 or (50 1 the maximal value 2^"^ is attained. Due to 


( |A2| ) we need to estimate the sum 

S„. = 


1 


2n+(n-l)/2 


E 






(A7) 
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and show that it can be equal to one by absolute value. First, let us consider the case of an odd n. From (501 
we have 






Furthermore, from (45 1 it is 


Sn = - 

2^i 


n(-r+-r)-n(4 

_k=l k=l 


10 ( 1 ) _ ■ * 0 ( 2 ) 


Taking into account that each term in these products can be written as 




and using the relations (A4 1 , Sn can be simplified to 

/ n n \ 

S'**=T (n ^ n 

* \fc=l k=l J 

where = 6k — 7^14:, (p'/^ = pk + /4. It is clear that the imaginary part of the sum Sn takes its maximal 
absolute value 1 when, for example, 9k = Pk = 0 for k = 1,... ,n. 

Now we consider the case of an even n. The coefficients in this case come from ([47]) 






and the sum SniAl \ becomes 


‘^®® ~ ^ ^®® ' llfc=l I lifc=l 


< ) re=i K + € ’ ) - WU 


er - e 


According to (A4 1 Sn can be simplified to 


Q _ 

- V2 


* nLi cos T riLi sin ) cos On + (JlLi cos d'l ± i OLi sin 9^ ) sin 0, 


jn-l n! 


in—1 


\n-l qH 


where 9'f. = 9k — vr/d, 9'1 = + vr/d and p = pi + ... Pn + {n — l)7r/4. The imaginary part of Sn is 

(when (/? = 0) 


|Im(Sn)| cos 0^=1 cos(0fc - 7r/4) ± sin 0^=1 sinC^fc + vr/d) 


It is clear that this expression takes its maximal value 1 when 9k = it( 4, k = 1,..., n — 1, and 9n = ±7r/4. 
This completes the proof. ■ 
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